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ORBITAL  TRANSFER  IN  MINIMUM  TIME 
V.  E.  BLEICK,*  U.  S.  Naval  Postgraduate  School 

1 .  Introduction.  The  problem  of  orbital  transfer  discussed 
here  is  that  of  scheduling  the  direction  p  of  constant  momentum 
thrust  of  a  rocket,  which  loses  mass  at  a  constant  rate,  so  that 
it  transfers  to  an  earth  satellite  orbit,  with  known  elements  of 
time,  position  and  velocity,  in  a  minimum  time  T  after  launching 
of  the  rock'et.  The  launching  conditions  are  assumed  to  be  fixed. 
This  situation  is  illustrated  in  Figure  1  for  the  case  of  a  cir¬ 
cular  orbit.  The  sector  angle  B  at  which  the  rocket  enters  orbit 
will  be  called  the  rendezvous  angle.  To  aid  the  discussion  imagin¬ 
ary  physical  rendezvous  of  the  rocket  and  satellite  is  assumed  to 
occur  at  this  angle.  The  time  of  rocket  launch  to  achieve  actual 
physical  rendezvous  can  be  determined,  of  course,  only  after 
both  of  the  unknowns  T  and  B  have  been  found.  The  problem  is  set 
up  as  a  calculus  of  variations  problem  of  the  Lagrange  type,  and 
is  solved  by  an  iterative  process  in  which  an  initial  approxima¬ 
tion  to  the  angle  B  is  estimated. 

A  non-rotating  Oxy  rectangular  coordinate  system  with  origin 
at  the  earth's  center  is  used.  The  coordinates  and  velocity  com¬ 
ponents  of  the  rocket  and  target  satellite  are  denoted  by  x,y,u,v 
and  X,Y,U,V  respectively.  For  simplicity  the  equations  of  motion 
of  rocket  and  target  will  be  written  in  a  "non-dimensional"  form 
by  the  use  of  suitable  units.  The  unit  of  length  is  taken  as  the 
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Fig.  1 .  Orbital  transfer. 


earth's  equatorial  radius,  =  20,925,000  feet.  The  unit  taken 
for  time  t  is  the  time  required  by  a  hypothetical  eart^i  satel¬ 
lite,  in  equatorial,  circular,  vacuum,  sea  level  orbit,  to 
traverse  a  sector  of  one  radian.  This  unit  of  time  is  YR~7g  = 
13.459  minutes,  where  g  =  32.086  ft. /sec.  is  the  acceleration 
of  gravity  at  the  equator.  The  unit  of  velocity  is  then  the  speed 
of  this  hypothetical  satellite.  These  units  of  length  and  time 
will  always  be  understood,  unless  other,  more  conventional  units 
are  specifically  mentioned. 

2.  Statement  of  the  problem.  The  equations  of  motion  of  the 
rocket,  in  terms  of  the  specified  units  of  length  and  time,  are 

(1)  ■  cp^  =  u  -  g^  -  a  cos  p  =  0 

cp2=v-g2-a  sinp  =  0 
cp^  =  X  -  u  =  0 
<^54  =  y  -  V  =  0 

where  g^  =  -x/r  ,  g2  =  -y/r  ,  r  =  x  +  y  ,  a  =  cm/g(l  -  mt) , 
where  m  is  the  constant  fraction  of  initial  gross  rocket  mass 
lost  per  unit  of  time,  c  is  the  constant  speed  of  the  emitted 
rocket  gases,  and  g  is  the  acceleration  of  gravity  at  the  equator. 
The  fixed  initial  conditions  of  the  rocket  trajectory  are  taken  as 

(2)  x(0)=0,  u(0)=V^='VqCOs6 

y(0)  =  1  ,  v(0)  =  V2  =  Vq  sin  6  . 
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The  terminal  point  of  the  rocket  trajectory  is  variable  with 


(3)  x(T)  =X(T),  y(T)  =  Y(T) ,  u(T)  =  U(T) ,  v(T)  =  V(T) . 


It  is  also  assumed  that  the  rocket  thrust  is  turned  off  abruptly 
at  time  T.  This  discontinuity  vill  lead  to  a  trivial  steering 
corner  in  the  calculus  of  variations  problem.  Note  that  (l)  and 
(3)  imply  that  U(T)  =  -[X/R^]^  =  (T)  and  V(T)  =  -[i/R^]^  =  g2(T) 

where  R  =  |T  X  +  ^  Y| . 

The  problem  is  to  choose  the  control  variable  p  to  effect 


r 

Jo 


orbital  transfer  with  L  dt  minimized,  and  to  determine  the 
corresponding  rocket  trajectory.  This  problem  is  equivalent  to 
the  Lagrange  calculus  of  variations  problem  of  requiring  the 
integral 


(4) 


tT 

""  Jo 


(1  +  X.cp^  +  HCP2  +  ncp^  +  pcp^)dt 


to  be  stationary.  In  (4)  the  equations  (l)  are  regarded  as  con¬ 
straints  with  X(t),  u(t)  ,  rr(t)  ,  p(t)  introduced  as  continua  of 
Lagrangian  multipliers  [l ] ,  If  the  time  coordinate  of  the  varied 
terminal  point  is  taken  as  T  +  AT ,  the  vanishing  first  variation^ 
[2]  of  I  is 

(5)  6I  =  j:  [x(6u  -  -  g-jy^y  +  a  sin  p6p)  +  tt(6x  -  6u) 

+  n(6v  -  g2x^x  -  S2y^y  ~  P^P)  +  P(6y  -  6v)  ]dt 

f  T+AT 

+  /m  [1  +  ^6u  +  n6v  -  (X.  cos  p  +  pi  sin  p)6a]dt  =  0 
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vhere  the  finite  variation  6a,  6u,  6v  terms  in  the  integral  from 
T  to  T  +  AT,  created  by  thrust  termination  at  time  T  on  the  un¬ 
varied  trajectory  and  at  T  +  AT  on  the  varied  trajectory,  cancel 
On  integrating  by  parts  one  obtains 

fT  , 

(6)  6l  =  [  X6u  +  ia6v  +  tt6x  +  p6y]^  -  [(X  +  tt)6u  +  (la  +  p)6v 

+  (p  +  g-|y^  +  g2y^^)^y 
+  a(|-i  cos  p  -  X  sin  p)6p]dt  +  AT  =  0. 

The  variations  of  the  dependent  coordinates  at  the  variable 

terminal  point  must  be  taken  [.2]  as 

(7)  6u(T)  =  (U  -  u)^AT  =  -(a  cos  p)^AT,  6x(T)  =  (X  -  x)^AT  =  0 

6v(T)  =  (V  -  v)^AT  =  -(a  sin  p),j,AT,  6y(T)  =  (Y  -  y)^AT  =  0, 

Substitution  of  (7)  into  (6),  and  application  of  the  fundamental 
lemma  of  the  calculus  of  variations  and  the  theory  of  ordinary 
extrema,  gives  the  Euler  equations 

(8)  X  +  TT  =  0 

+  p  =  0 

"  +  =  0 

P  +  g-t/  +  gj/  =  ® 
tan  p  =  n/X 

and  the  transversality  condition 

[a(X  cos  p  +  H  sin  p)  Jrj,  =  1. 
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The  first  four  homogeneous  equations  of  the  Euler  equations  (8) 
constitute  the  adjoint  system  [3]  of  the  system  of  variation 
equations 

(10)  6cp^  =  6cp2  =  ^<¥'3  =  ^^4  =  0 

which  are  the  coefficients  of  X,  tt,  p  in  (5)  equated  to  zero. 

The  adjoint  system  has  a  matrix  of  coefficients  which  is  the  neg¬ 
ative  transpose  of  that  of  (10).  The  last  of  the  Euler  equations 
(8)  requires  that  the  control  variable  p  be  adjusted  so  that 
a  =  a(T  cos  p  +  T  sin  p) ,  which  is  proportional  to  the  rocket 
thrust,  is  continually  parallel  to  the  adjoint  vector  A  =  i  X  +  3  p. 
The  transversality  condition  (9) ,  which  may  be  written 
(a*A)^  =  1  >  0,  requires  that  a  and  A  have  the  same  sense.  Since 
it  is  only  the  ratio  of  H  to  X  which  determines  p,  it  is  a  trivial 
matter  to  scale  them  to  satisfy  the  magnitude  requirement  of  (9). 

A  solution  of  the  problem  obtained  from  (I)  and  (8)  guarantees 
a  stationary  time  of  transfer.  The  nature  of  the  problem  is  such 
that  this  stationary  time  is  a  minimum  time. 

3.  Numerical  solution.  There  is  a  constructive  aspect  of  a 
modification  of  equation  (6),  first  used  by  Bliss  [4]  in  his  work'* 
on  differential  corrections  in  ballistics,  and  applied  recently 
by  Faulkner  [5]  in  an  iterative  fashion  in  optimum  control  prob¬ 
lems.  To  find  the  desired  modification  of  (6),  assume  that  a 
solution  of  the  systems  (I)  and  (8)  has  been  obtained,  which  does 
not  necessarily  satisfy  the  terminal  conditions  (3).  Using  this 
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solution  and  holding  T  constant,  consider  the  variation  of  the 
vanishing  integral 


(11) 


+  +  TTCp^  +  pCP4)dt  =  0 


with  the  terminal  constraints  (3)  removed,  so  that  the  terminal 
variations  6u(T) ,  6v(T) ,  6x(T)  and  6(T)  become  free.  Since  X,  n, 
TT,  p  satisfy  the  adjoint  system,  and  since  there  is  now  no  steer¬ 
ing  corner  due  to  thrust  termination,  one  obtains 


(12)  fx6u  +  u6v  +  tt6x  + 


r  T 

p^yjf  =  Jo  sin  p  +  U  cos  p)6pdt 

TT  ^ 

=  Jq  a* ( da/dp) 6pdt 


where  a  =  a(T  cos  p  +  J  sin  p) .  Equation  (12),  which  is  the 
desired  modification  of  (6),  is  called  the  fundamental  formula 
by  Bliss  [4],  but  is  also  known  under  the  generic  name^ of  Green's 
formula  [3].  By  the  use  of  (12)  it  is  possible  to  generate  the 
control  parameters  of  a  varied  trajectory  which,  hopefully, 
comes  closer  to  satisfying  the  desired  terminal  conditions  (3). 

To  do  this,  assume  that  the  adjoint  system  has  been  solved  to 
obtain  a  fundamental  set  of  four  linearly  independent  solutions 
given  by  the  rows  of 

(13)  Transpose  of  B(t)  =  [X^(t)  U^(t)  TT^(t)  P^(t)]  i  =  1,2, 3, 4 

where  B(0)  =  I  is  the  identity  matrix.  The  solution  A  =  T  X  +  J  pi 
of  the  adjoint  system,  required  to  satisfy  the  last  of  the  Euler 
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equations  (8) ,  is  taken  as  the  linear  combination 


(14)  \  ^^4 

li  =  +  1^2  +  naU^  + 

so  that  the  control  angle  p  is  determined  by 

(15)  tan  p  =  (u^  +  l^i2  +  +  np^)/(X^  +  1^2  +  m^3  +  nX^) 

and  its  variation  by 

(16)  6p  =  [(Xu2-uX2)61  +  (Xn^-uX^)6m  +  ( X^^-^iX^) 6n]/( X^  +  . 

When  (13)  and  (I6)  are  substituted  into  (12)  there  results  the 
system  of  Green's  formulae 

(17)  [6u  6v  6x  6y]^B(T)  =  [O  6l  6m  6n]A 
where  the  elements  of  the  matrix  A  are 

(18)  a.^  =  /o^a(Xu.  -  uX. ) (Xu^  -  ^X^)dt/(x2  + 

The  coordinates  of  the  terminal  point  of  the  varied  trajectory 
at  time  T  +  AT  may  be  taken  as  [u+Au  v+Av  x+Ax  y+Ay],p  where 

(19)  [Au  Av  Ax  Ay]^  =  [6u  6v  6x  6y]^  +  [u  v  x  y],pAT. 

In  an  effort  to  make  this  new  terminal  point  come  closer  to  satis¬ 
fying  the  terminal  conditions  (3),  one  may  take 

(20)  [Au  Av  Ax  Ay]^  =  [U-u  V-v  X-x  Y-y]^  +  [U  V  X  Y]^AT. 
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Substitute  (19)  and  (20)  into  (1?)  to  obtain 

(21)  [U-u  V-v  X-x  Y-y]^AT  +  [O  6l  6ra  6n]AB"^  (T)  =  [U-u  V-v  X-x  Y-y]q, 

as  the  system  of  equations  for  the  determination  of  AT,  6l,  6m,  6n 
on  the  varied  trajectory.  The  Faulkner  [5]  scheme  for  the  numerical 
solution  of  optimum  control  problems  may  nov  be  stated:  Make  an  ^ 
initial  guess  for  the  values  of  T,  1,  m,  n;  carry  out  a  simultan¬ 
eous  numerical  integration  of  the  systems  (1),  (8)  and  (18)  using 

the  control  variable  of  (15);  solve  the  system  (21)  for  AT  and 
the  changes  in  the  control  parameters;  iterate  until  convergence 
is  obtained.  This  program  may  be  carried  out  in  a  matter  of  seconds 
on  modern  digital  computers. 

To  give  an  example  of  this  control  optimization  in  the  present 
problem  a  circular  satellite  orbit  was  assumed  of  radius  R  = 

|i  X  +  ^  T|  =  1.075699,  corresponding  to  an  altitude  of  300  statute 
miles  above  sea  level.  For  this  orbit  |T  U  +  ^  V|  =  1 /YR  and 
I  i  U  +  j  V|  =  1/R^.  The  assumed  rocket  launching  velocity  was 
Vq  =  |T  ^  V2 I  =  0.585402,  launch  angle  6  =  0.928084  and  ren¬ 

dezvous  sector  B  =  0.153840  as  in  Figure  1.  Also  assumed  were 
c  =  10000.9  ft. /sec.  and  lii  =  0.00360583  sec.  ^  The  odd  appearance 
of  these  figures  is  related  to  the  difficulty,  explained  in  the 
next  section,  of  obtaining  the  initial  "guesses"  T  =  0.289725, 

1  =  -0.223125,  ra  =  -29.9875  and  n  =  19.0847  .  Vith  this  input  the 


process  converged  in  four  iterations  to  the  seven  significant 
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Pig.  2.  Trajectory  and  thrust  directions. 


figure  results  T  =  0.2894592,  1  =  0.1840054,  m  =  -108.94383, 
n  =  67.95886  and  B  =  0.1536015  .  Figure  2  shows  the  resulting 
trajectory  and  rocket  thrust  directions  for  equal  time  intervals, 
excepting  the  interval  terminating  in  transfer,  which  is  reduced 
by  one  half. 

4.  The  initial  guesses.  The  domain  of  convergence  of  the 
iteration  scheme  of  the  last  section  appears  to  be  rather  limited 
in  the  present  control  problem,  requiring  initial  guesses  for 
T,  1,  m,  n  which  make  [U-u  V-v  X-x  Y-y]rp  small  in  (21  )  .  In 
the  numerical  example  given  here,  where  the  angle  B  is  small  and 
radius  R  nearly  unity,  the  equations  (l)  can  be  linearized  and  an 

exact  solution  of  the  linearized  transfer  problem  used  to  supply 

> 

the  input  guesses  required  to  solve  the  non-linear  problem.  The 
gravitational  terras  g^  and  g2  of  (l)  were  replaced  by  the  linear 
terras  in  their  'faylor  expansions  at  x  =  0,  y  =  1.  The  systems  (1) 
and  (8)  then  become  uncoupled.  Their  solutions,  using  (2),  are 


(22) 


cos  p  sin  (t-w)dw  +  sin  t 
cos  p  cos  (t-w)dw  +  cos  t 

sin  p  sinh  Y^(t-w)dw  +  sink  Y^t  -  (l/)r2)cosh  Y^t  +  3/Y^ 

sin  p  cosh  Y^(t-w)dw  +  V2  cosh  Y^t  -  (l/Y2)sinh  Tst 


tan  p  =  [1  cosh  Y^t  -  (n/Y^)sinh  Y^t]/(cos  t  -  m  sin  t) . 
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It  will  be  of  no  avail  to  attempt  to  solve  (22)  and  (3)  for 
T,  1,  m,  n  by  the  Nevton-Raphson  method,  since  the  Newton-Raphson 
equations  are  in  fact  the  system  (21)  of  the  poorly  convergent 
iterative  routine  of  the  last  section.  A  substitute  procedure 
of  solving  (22)  and  (3)  for  ,  V2 ,  m  and  B  was  used.  Vhen  the 
transfer  problem  has  been  solved  for  a  sufficient  number  of  such 
sets  (V^ ,  V2 ,  m,  B) ,  a  basis  will  be  at  hand  for  obtaining  desired 
rocket  launching  conditions  by  interpolation  or  extrapolation.. 

It  can  now  be  revealed  that  the  numerical  example  of  the  last 

section  really  had  its  genesis  in  the  assumptions  B  =  n/l6,  9  =  tt/4, 

»  -1 

Vq  =  0.65,  m  =  0.0025  sec.  and  c  =  10000  ft. /sec.  It  was 
estimated  that  T  =  0.3,  1  =1.25,  m  =  -1.0  and  n  =  9.5  would 
satisfy  (22)  and  (3).  Using  a  cluster  of  five  closely  spaced 
points  around  (T,l,m,n),  a  single  application  of  regula  falsi  [6] 
was  made  in  the  hope  of  improving  these  values  of  T,  1,  m,  n.  The 
particular  cluster  chosen  gave  the  output  values  T  =  0.289725, 

1  =  -0.223125,  m  =  -29.9875,  n  =  19.0847,  which  will  be  recognized 
as  the  initial  guesses  of  the  last  section.  The  output^ residuals  were 
[X-x,  T-y,  U-u,  V-v]^  =  [0.041389,  -0.020160,  0.212093,  -0.221912]. 

These  residuals  were  then  processed  by  the  linear  system 

fT 

(23)  (X-x)^  =  AV^  sin  T  -  T(T)dB  +  q  a  cos  p  sin  (T-w)dw/c 

r  ^ 

(U-u)m  =  AV.  cos  T  -  V(T)dB  +  Ac /„  a  cos  p  cos  (T-w)dw/c 

r-  ''  /“T 

Y2(Y-y)^  =  AV2sinh  Y^T  +  Y^X(T)dB  +  ^cJq  a  sin  p  sinh  Y^(T-w)dw/c 
(V-v)^  =  AV2COsh  YIt  +  U(T)dB  +  AcJ'^'^  a  sin  p  cosh  Y^(T-w)dw/c, 
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derived  from  (22)  and  (3)  for  the  case  of  a  circular  orbit,  to 
obtain  c  +  Ac,  launching  conditions  +  AV^  and  V2  +  ^"^2' 
dezvous  sector  B  +  dB  which  will  give  small  residuals.  The  value 
of  c  +  Ac  turned  out  to  be  far  from  the  desired  10000  ft. /sec. 
Instead  of  changing  the  value  of  c,  the  transformation 

(24)  10000  In  [1  -  (m  +  Am)T]  =  (c  +  Ac)  In  (1  -  mT) , 

/■T 

which  leaves  the  integral  Jq  adt  invariant,  was  used  to  change  m 
and  cause  c  +  Ac  to  approach  10000  ft. /sec.  on  iterating  (24), 

^  (22)  and  (23).  Seven  iterations  produced  the  input  values  of  V.j  , 
,  c,  m  and  B  used  in  the  last  section. 

5.  The  problem  of  two  fixed  end  points.  The  problem  in  which 
the  conditions  at  both  ends  of  the  rocket  trajectory  are  fixed, 
or  the  problem  in  which  rocket  and  target  satellite  achieve 
actual  physical  rendezvous  at  a  specified  angle  B,  can  also  be 
solved  by  the  methods  presented  here,  but  with  greater  converg¬ 
ence  difficulties. 

Acknowledgement.  1  am  indebted  to  my  colleague.  Prof.  Frank 
D.  Faulkner,  for  suggesting  the  problem  of  this  paper  and  for 
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CDC-1604  ASSEMBLY  LANGUAGE  COMPUTER  PROGRAM 
Set  content  of  K1 00+1  «  SLJ  0  L3440  ZRO  0  00000 


KO 


KIO 


K20 


K30 


REM 

REM 

REM 

REM 

REM 

SIL 

1 

1 

NT 

COMPUTE  LINEARIZED  TRAJECTORY 
(81)  =  NUMBER  OF  TRAJECTORIES 
(B4)  =  NUMBER  OF  FALSI  ITERAT 
(A)  =  FIRST  DEPARTURE  ANGLE 
(Q)  =  DELTA  DEPARTURE  ANGLE 
STORE  NO.  OF  TRAJS. 

STA 

0 

THETATEM 

STORE  1ST  DEPART  ANGLE 

STQ 

0 

DELTHETA 

STORE  DELTA  DEPART  ANGLE 

LDA 

0 

UNITY 

FDV 

0 

R 

RECIP.  ORBIT  RADIUS 

SLJ 

4 

SQROOT 

RET. JUMP  TO  SQUARE  ROOT 

STA 

0 

V 

STORE  TARGET  SPEED 

FMU 

0 

V 

FDV 

0 

R 

STA 

0 

VSQDR 

STORE  TARGET  ACCEL. 

LDA 

0 

B 

TARGET  SECTOR  AT  RENDEVOUS 

SLJ 

4 

TRIG+70 

RET. JUMP  TO  TRIGONOMETRIC  SR. 

STA 

0 

COSB 

FMU 

0 

R 

STA 

0 

CAPYFIN 

TARGET  COOR.  AT  RENDEVOUS 

LDA 

0 

B 

SIU 

4 

K20  +  6 

STORE  B4 

SLJ 

4 

TRIG 

STA 

0 

SIN8 

FMU 

0 

R 

STA 

0 

CAPXFIN 

TARGET  COOR.  AT  RENDEVOUS 

LDA 

0 

COSB 

FMU 

0 

V 

STA 

0 

CAPUFIN 

TARGET  VEL.COMP.  AT  RENDEVOUS 

LAC 

0 

SIN8 

FMU 

0 

V 

STA 

0 

CAPVFIN 

TARGET  VEL.COMP.  AT  RENDEVOUS 

LAC 

0 

COSB 

FMU 

0 

VSQDR 

STA 

0 

CAPVDFN 

TARGET  ACC. COMP.  AT  RENDEVOUS 

LAC 

0 

SINB 

FMU 

0 

VSQDR 

STA 

0 

CAPUDFN 

TARGET  ACC. COMP.  AT  RENDEVOUS 

LDA 

0 

REARTH 

EARTH  RADIUS  AT  EQUATOR 

FDV 

0 

GACCEL 

ACCEL .GRAVITY  AT  EQUATOR 

SLJ 

4 

SQROOT 

RECIPROCAL  SHULER  FREQ. 

FMU 

0 

MOOT 

STA 

0 

OMEGA 

NON-DIMEN.  MASS  LOSS  RATE 

LDA 

0 

C 

ROCKET  NOZZLE  VELOCITY 

FMU 

0 

MOOT 

C  MOOT 

FDV 

0 

GACCEL 

STA 

0 

A 

NON-DIMEN.  ACCELERATION 

LDA 

0 

THETATEM 

BEGIN  OUTER  LOOP 

STA 

0 

THETA 

SET  DEPARTURE  ANGLE 

ENI 

4 

0 

SET  NO.  OF  FALSI  ITERATES 

SLJ 

4 

TRIG+70 

COSINE  THETA 

FMU 

0 

VSTART 

STA 

0 

V  1 

INITIAL  HORIZONTAL  VELOC. 

LDA 

0 

THETA 

SLJ 

4 

TRIG 

SINE  THETA 

FMU 

0 

VSTART 

STA 

0 

V2 

INITIAL  VERTICAL  VELOC. 

ENI 

3 

23 

ENI 

0 

0 

LDA 

3 

EL1INIT 

SET  FALSI  COOR.  QUINTUPLE 

STA 

3 

ELI 

UP 

3 

7-1 

SLJ 

4 

70007 

PRINT  COORS.  DECIMAL.  B3=0 

STA 

0 

ELIINIT 

PARAMETER  WORD 

LESS  QNE 


13 


K40 


K50 


K60 


K70 


KlOO 


ZRO 

0 

T5INIT 

ENI 

b 

0 

SET  B5=0  FOR  K40+3 

SLJ 

4 

700C7 

PRINT  COORDINATES  OCTAL 

ZRO 

0 

EL  UNIT 

PARAMETER  WORD 

ZRO 

0 

T5INIT 

SLJ 

4 

XYUV 

SET  FALSI  FUNC.  QUINTUPLE 

ZRO 

0 

0 

SLJ 

4 

70007 

PRINT  FUNCTIONS  DECIMAL 

ZRO 

0 

0 

STA 

0 

DELXF IN 

PARAMETER  WORD 

ZRO 

0 

DELVFIN 

LOA 

0 

DELXFIN 

STA 

5 

OELXl 

LDA 

0 

DELYFIN 

STA 

5 

OELYl 

LDA 

0 

OELUFIN 

STA 

5 

0ELU1 

LDA 

0 

DELVF  IN 

STA 

5 

DELVl 

INI 

5 

5 

INI 

3 

3 

I  SK 

3 

23 

SLJ 

0 

K40 

ENI 

0 

0 

BEGIN  INNER  LOOP  OR  SLJ.O  KlOl  01 

SLJ 

4 

MATRIX 

INVERT  DELXYUV  MATRIX 

SLS 

0 

K50+2 

SINGULAR  MATRIX  ALARM  HALT 

ZRO 

0 

0 

ZRO 

0 

0 

ZRO 

0 

5 

CODE  TO  INVERT  =  S 

ZRO 

0 

5 

L  =  S  ROWS  IN  DELXYUV  MATRIX 

ZRO 

0 

DELXl-1 

ADD.  OF  DELXYUV  MATRIX  =  DELXl-1 

ZRO 

0 

0 

N=0.  WORD  NOT  USED 

ZRO 

0 

0 

ZRO 

0 

0 

M=0 

ZRO 

0 

DELXl+174 

ADO.  OF  INVERSE  =  DELXl+174 

ENI 

0 

0 

SLJ 

4 

MATRIX 

INVERSE  TIMES  LMNT  MATRIX 

ZRO 

0 

K60 

WORD  NOT  USED 

ZRO 

0 

0 

ZRO 

0 

0 

ZRO 

0 

4 

MATRIX  MULT.  CODE  =  4 

ZRO 

0 

5 

L  =  S  ROWS  IN  INVERSE 

ZRO 

0 

DELXl+174 

ADD.  OF  INVERSE  =  DELXl+174 

ZRO 

0 

5 

N=S  COLUMNS  IN  INVERSE 

ZRO 

0 

ELI 

ADDRESS  OF  LMNT  MATRIX  =  ELI 

ZRO 

0 

4 

M=4  COLS.  IN  MATRIX  PRODUCT 

ZRO 

0 

ELl+24 

ADDRESS  OF  PROD.  =  ELl+24 

LDA 

3 

EL2 

SHIFT  COORS.  AND  FUNCS.  B3=0 

STA 

3 

ELI 

LDA 

3 

DELX2-1 

STA 

3 

DELXl- 1 

ISK 

3 

23 

SLJ 

0 

/-2 

ENI 

3 

20 

ENI 

5 

3 

FOR  USE  IN  K70+2 

SLJ 

4 

XYUV 

COMPUTE  NEW  DELXYUV  FUNCS. 

ZRO 

0 

0 

LOA 

5 

DELXFIN 

STORE  NEW  FUNCS. 

STA 

5 

0ELX5 

UP 

5 

/-I 

. 

SLJ 

4 

70007 

PRINT  NEW  COORS.  DEC.  WITH  PANEL 

SAU 

0 

ELS 

PARAMETER  WORD 

ZRO 

0 

T5 

EXE 

0 

70 

CLEAR  ARITHMETIC  ERRORS 

SLJ 

4 

70007 

PRINT  NEW  COORS.  OCTAL 

ZRO 

0 

ELS 

PARAMETER  WORD 

ZRO 

0 

TS 

ENI 

3 

0 

CLEAR  83 

SLJ 

4 

70007 

PRINT  NEW  FUNCS.  DECIMAL 

STA 

0 

DELXFIN 

PARAMETER  WORD 

14 


ZRO 

0 

DELVFIN 

UP 

4 

K50  +  1 

LDO 

0 

ELS 

LDA 

0 

UNITY 

SLJ 

4 

POLAR+130 

STQ 

0 

OQ 

ENI 

T 

2 

LDA 

1 

ELS 

STA 

1 

EL 

UP 

1 

/-I 

LAC 

0 

V 

FMU 

0 

TAUFIN 

FDV 

0 

R 

FAD 

0 

B 

STA 

0 

S 

Kl  10 

SLJ 

4 

70007 

ZRO 

0 

0 

SAU 

0 

CAPXFIN 

ZRO 

0 

LAM 

SLJ 

4 

70007 

ZRO 

0 

0 

STA 

0 

NT 

ZRO 

0 

S 

SLJ 

4 

70007 

ZRO 

0 

0 

ZRO 

0 

NT 

ZRO 

0 

DELTHETA 

LDA 

0 

THETA 

FAD 

0 

DELTHETA 

STA 

0 

THETATEM 

RSO 

0 

NT 

K120 

AJP 

1 

K20+5 

SLS 

0 

BEGIN 

XYUV 

SLJ 

0 

0 

LDA 

3 

T1 

STA 

0 

TAUFIN 

LAC 

0 

OMEGA 

FMU 

0 

TAUFIN 

FAD 

0 

UNITY 

STA 

0 

STORl 

FDV 

0 

OMEGA 

STA 

0 

STOR2 

FMU 

0 

STOR2 

FMU 

0 

STOR2 

FMU 

0 

STOR2 

FMU 

0 

STOR2 

FMU 

0 

ERRCOF 

K130 

FDV 

0 

TAUFIN 

STA 

0 

STOR3 

LDA 

0 

STORl 

SLJ 

4 

LOG+66 

FMU 

0 

ST0R3 

SLJ 

4 

SQROOT 

SLJ 

4 

SQROOT 

ZRO 

0 

0 

STA 

0 

DELTAU 

SLJ 

4 

EQUAL 

STA 

0 

DELTl 

ENA 

0 

4 

STA 

0 

RUNKTl 

ENI 

1 

14 

ENA 

0 

0 

STA 

1 

TAUl 

K140 

UP 

1 

/-I 

EN! 

0 

0 

LDA 

0 

TAUFIN 

SLJ 

4 

TRIG 

S  *  A 

0 

SINT 

FMU 

0 

VI 

FSB 

0 

CAPXFIN 

END  INNER  LOOP  (  OR  SLJ  0 


DEPARTURE  THRUST  ANGLE 
MOVE  ELS,  EM5  AND  EN5 


(-VT) 

(-VT/R) 

B-( VT/R) 

S  =  INITIAL  TARGET  SECTOR 
ALARM  REENTRY  FROM  LOG  S.R. 
DECIMAL  DUMP  WITH  PANEL 
PARAMETER  WORD 

DECIMAL  DUMP 

PARAMETER  WORD 

OCTAL  DUMP 

PARAMETER  WORD 


SET  NEW  THETA 

END  OUTER  LOOP.  (SEE  K440) 
HALT.  START  RENDEVOUS  PROG. 
ENTER  DELXYUV  SUBROUTINE 


(STORl )  =  1  -  WT 

(ST0R2)  =  (1-WT)/W 

SQUARE 

CUBE 

FOURTH 

FIFTH 

MULT.  BY  -15.B-12 


1-WT 

LOG(l-WT)  TO  BASE  E 


COMPUTED  TAU  INCREMENT 
EQUALIZE  DELTAU  AND  SET  B2 
SET  EQUALIZED  DELT1 

SET  GILL  FOR  4-  INTEGRALS 


CLEAR  RUNKTl  AREA 


SIN  T 

VI  SINT 
VI  SINT  -  X 


-  15  - 


L3440  ) 


K150 


KI60 


K170 


K200 


STA 

0 

DELXFIN 

LOA 

0 

TAUFIN 

SLJ 

4 

TRI  G+70 

COS  T 

STA 

0 

COST 

FMU 

0 

V  1 

VI  COST 

FSB 

0 

CAPUFIN 

VI  COST  -  U 

STA 

0 

DELUFIN 

LDA 

0 

TAUFIN 

FMU 

0 

ROOT2 

T  TIMES  SQ.  ROOT  2 

SLJ 

4' 

EXP+107 

ZRO 

O' 

0 

STA 

0 

COSHTR2 

STORE  HYP.  COSINE 

FMU 

0 

ROOTHALF 

COSHTR2/SQ.ROOT  2 

STA 

0 

STORl 

LDA 

0 

TAUFIN 

FMU 

0 

R00T2 

SLJ 

4 

EXP+71 

STA 

0 

SINHTR2 

STORE  HYP.  SINE 

FMU 

0 

V2 

V2  SINHTR2 

FSB 

0 

STORl 

FAD 

0 

THREOR2 

ADD  (3/SQ.ROOT  2) 

STA 

0 

STOR  1 

LAC 

0 

CAPYFIN 

(-Y) 

FMU 

0 

ROOT2 

FAD 

0 

STORl 

STA 

0 

DELYFIN 

LDA 

0 

COSHTR2 

FMU 

0 

V2 

STA 

0 

STORl 

(STORl)  =  V2  COSHTR2 

LAC 

0 

SINHTR2 

(-SINHTR2) 

FMU 

0 

ROOTHALF 

FAD 

0 

STORl 

FSB 

0 

CAPVFIN 

STA 

0 

DELVFIN 

LDA 

3 

EL  1 

FMU 

3 

ELI 

ELI  SQUARED 

FAD 

0 

UNITY 

1  +  ELI  SQ. 

SLJ 

4 

SQRCOT 

ZRO 

0 

0 

STA 

0 

HYP 

LDA 

0 

A 

FOV 

0 

HYP 

STA 

0 

STORl 

(STORl)  =  A/HYP  =  A  COS  P 

FMU 

0 

SINT 

STA 

0 

SINDOT 

A  COSP  SINT 

LDA 

0 

STORl 

A  COSP 

FMU 

0 

COST 

STA 

0 

COSDOT 

A  COSP  COST 

LDA 

0 

STORl 

A/HYP 

FMU 

3 

ELI 

A  SINP 

FMU 

0 

SINHTR2 

STA 

0 

SINHDOT 

A  SINP  SINHTR2 

LDA 

0 

STORl 

A/HYP 

FMU 

3 

ELI 

A  SINP 

FMU 

0 

COSHTR2 

STA 

0 

COSHDOT 

A  SINP  COSHTR2 

SLJ 

4 

ALPHA 

SET  UP  GILL  ROUTINE 

ZRO 

0 

RUNKTl 

PARAMETER  WORD 

ZRO 

0 

DERIV1 

UP 

2 

/+! 

JUMP  IF  B2  NOT  =  ZERO 

SLJ 

0 

/+3 

62=0.  END  OF  INTEGRATION 

6NI 

0 

0 

62  NOT  =  ZERO 

SLJ 

4 

ALPHA+ 1 

INTEGRATE  AGAIN 

ENI 

0 

0 

SLJ 

0 

/-2 

JUMP  TO  TEST  B2 

LDA 

0 

DELXFIN 

FAD 

0 

SIN 

STA 

0 

DELXFIN 

LDA 

0 

DELUFIN 

FAD 

0 

COS 

-  16  - 


K2t0 


DERIV1 


K220 


K230 


K240 


K250 


STA 

LDA 

FAD 

FDV 

STA 

LDA 

FAD 

STA 

SLJ 

LAC 

FHU 

FAD 

STA 

LDA 

FDV 

STA 

LDA 

FSB 

STA 

ENI 

SLJ 

STA 

LDA 

SLJ 

ZRO 

STA 

LDA 

FMU 

SLJ 

STA 

LDA 

FMU 

SLJ 

STA 

LDA 

SLJ 

ZRO 

FMU 

STA 

LDA 

SLJ 

FSB 

STA 

FMU 

STA 

LDA 

FMU 

SLJ 

ZRO 

FMU 

FDV 

STA 

LDA 

FMU 

SLJ 

FMU 

FS8 

STA 

FMU 

FAD 

SLJ 

STA 

LDA 

FDV 

STA 

LDA 

FDV 

STA 

FMU 

FMU 


0 

DELUFIN 

0 

OELYFIM 

0 

SINH 

0 

ROOT2 

0 

DELYFIN 

0 

OELVFIN 

0 

COSH 

0 

OELVFIN 

0 

XYUV 

TO  EXIT  OF  SUBROUTINE 

0 

TAUl 

0 

OMEGA 

(-WT) 

0 

UNITY 

0 

UNMINWT 

1-WT 

0 

A 

0 

UNMINWT 

0 

ADOMINWT 

A/(  1-WT) 

0 

TAUFIN 

0 

TAUl 

(STOR 1 )  =  TAUFIN  -  TAU 

0 

STORl 

0 

0 

4 

TRIG 

(SINT)  =  SINE  OF  (STOR 

0 

SINT 

0 

STORl 

4 

TRIG+70 

0 

0 

0 

COST 

0 

STORl 

0 

ROOT2 

4 

EXP+71 

0 

SINHTR2 

0 

STORl 

0 

R00T2 

4 

EXP+107 

0 

C0SHTR2 

0 

TAU  1 

4 

TRIG 

SIN  TAUl  =  -LAM3 

0 

0 

3 

EMI 

0 

ST0R2 

(ST0R2)  =  -  EMI  LAM3 

0 

TAUl 

4 

TRIG+70 

COS  TAUl  =  LAMl 

0 

ST0R2 

LAMl  +  EMI  LAM3 

0 

LAM 

LAM  =  LAMl  +  EMI  LAM3 

0 

LAM 

0 

ST0R3 

{STOR3)  =  LAM  SQUARED 

0 

TAUl 

0 

ROOT2 

4 

EXP+71 

SINH  TAU1R2 

0 

0 

3 

ENI 

0 

ROOT2 

(  -  ENI  MU4  ) 

0 

STOR2 

.  (STOR 2)  =  -  ENI  MU4 

0 

TAUl 

0 

ROOT2 

4 

EXP+107 

COSH  TAU1R2  =  MU2 

3 

ELI 

ELI  MU2 

0 

STOR2 

ELI  MU2  +  ENI  MU4  =  MU 

0 

MU 

0 

MU 

0 

STOR3 

LAM  SO.  +  MU  SQ. 

4 

SQROOT 

0 

HYP 

0 

MU 

0 

HYP 

0 

ST0R3 

(STOR3)  =  SINP 

0 

Q 

LAM 

HYP 

0 

.  ST0R4 

(STOR4)  »  COSP 

0 

SINT 

COSP  SINT 

0 

ADOMINWT 

-  17  - 


STA 

0 

SINDOT 

LDA 

0 

ST0R4 

COSP 

FMU 

0 

COST 

COSP  COST 

FMU 

Q 

ADOMINWT 

STA 

5 

COSDOT 

LDA 

0 

ST0R3 

SINP 

FMU 

0 

SINHTR2 

SINP  SINHTR2 

FMU 

0 

ADOMINWT 

STA 

0 

SINHDOT 

LDA 

0 

STOR3 

SINP 

FMU 

0 

COSHTR2 

SINP  COSHTR2 

FMU 

0 

ADOMINWT 

STA 

0 

COSHDOT 

EXIT  FROM  DERIVl 

K260 

SLJ 

0 

ALPHA+2 

ZRO 

0 

0 

EQUAL 

SLJ 

0 

0 

S.R.  TO  EQUALIZE  DELTAU 

LDA 

0 

TAUFIN 

FDV 

0 

DELTAU 

CONVERT  DELTAU  TO  FIXED  PT. 

SLJ 

4 

FIX  IT 

STA 

0 

INTSIGN 

STORE  INTEGRAL  PART 

STQ 

0 

FRACSIGN 

STORE  FRACTIONAL  PART 

AJP 

1 

/+5 

JUMP  IF  A  NOT  =  ZERO 

QJP 

0 

/+3 

A=0.  JUMP  IF  Q  =  0 

LDA 

0 

TAUFIN 

A=0.  Q  NOT  =  0 

STA 

0 

DELTFLSG 

6NI 

2 

1 

SLJ 

0 

/-H2 

ENA 

0 

0 

A  =  Q  =  0 

STA 

0 

DELTFLSG 

K270 

ENI 

2 

0 

SLJ 

0 

/+10 

AJP 

2 

/+1 

A  NOT  0.  JUMP  IF  A  GREATER  THAN  0 

SCM 

0 

MASK 

A  LESS  THAN  0.  ABSO  INTEGER  TO  A 

INA 

0 

1 

ABSO.  INTEGER  +  1  TO  A 

SAU 

0 

/+! 

ENI 

2 

0 

ABSO.  INTEGER  +  1  TO  B2 

SCA 

1 

2057 

LRS 

0 

57 

ENA 

1 

0 

LLS 

0 

44 

STA 

0 

FLOIPLON 

LDA 

0 

TAUFIN 

FDV 

0 

FLOIPLON 

STA 

0 

DELTFLSG 

ENI 

0 

0 

K300 

SIL 

2 

FIXIPLON 

USED  IN  FINAL  PROG. 

LDA 

0 

DELTFLSG 

SLJ 

0 

EQUAL 

TO  EXIT  OF  S.R. 

ZRO 

0 

0 

FIXIT 

BSS 

44 

Bl=5  B6=65302  P=460 

R 

DEC 

1.075698925  NON-DIMEN.  ORBIT  RADIUS 

V 

BSS 

1 

NON-OIMEN.  TARGET  SPEED 

VSQDR 

BSS 

1 

NON-DIMEN.  TARGET  ACCEL. 

B 

OCT 

1775622077325042  TARGET  SECTOR  AT  RENDEVOUS  = 

COSB 

BSS 

1 

SINB 

BSS 

1 

REARTH 

DEC 

20925000. 

EARTH  EQUATOR  RAD.  IN  FT. 

GACCEL 

DEC 

32.086 

EQUATOR  GRAVITY  IN  FT/SEC/SEC 

RMS 

BSS 

1 

OMEGA 

BSS 

1 

NON-DIMEN.  MASS  LOSS  RATE 

THETATEM 

BSS 

1 

VI 

BSS 

1 

NON-DIMEN.  INIT.  HORIZONTAL  VEL. 

V2 

BSS 

1 

NON-DIMEN.  INIT.  VERTICAL  VEL. 

ERRCOF 

DEC 

-15.B-12 

(-15/4096  DEC.)  USED  TO  FIND  DELTAU 

ROOTHALF 

OCT 

2000552023631500 

SINT 

BSS 

1 

COST 

BSS 

1 

SINHTR2 

BSS 

1 

C0SHTR2 

BSS 

1 

UNMINWT 

BSS 

1 

STORE  FOR  (l-WT) 

18 


EL  1 INIT 

EMI  INI T 

ENIINIT 

T1  INIT 

EL2INIT 

EM2INIT 

EN2INIT 

T2INIT 

EL3 INIT 

EM3INIT 

EN3INIT 

T3INIT 

EL4 INIT 

cMUINIT 

EN4INIT 

T4INIT 

EL5INIT 

EM5INIT 

EN5INIT 

T5INIT 

RUNKTl 

CELT  1 

TAU1 

SINOOT 

SIN 

COS  DOT 
COS 

SINHOOT 

SINH 

COSHDOT 

COSH 

K440 


EXP 

SOROOT 

TRIG 

AD^MINWT 

R00T2 

TWO 

THREE 

MASK 

TENGRANO 

COMPRMS 

BEGIN 


LIO 


BSS 

2 

OCT 

2001500400000000 

DEC. VALUE  = 

+1.251953125 

OCT 

5776377777777777 

DEC. VALUE  = 

-1.000000000 

OCT 

2004460000000000 

DEC. VALUE  = 

+9.500000000 

OCT 

1776463 146314631 

DEC. VALUE  = 

+0.3000000000 

OCT 

2001500300000010 

DEC. VALUE  = 

+  1  .25  1464844 

OCT 

5776377777777777 

DEC. VALUE  = 

-1.000000000 

OCT 

2004460100000000 

DEC. VALUE  = 

+9.503906250 

OCT 

1776461700000000 

DEC. VALUE  = 

+0.2987060547 

OCT 

200 1500200000000 

DEC. VALUE  = 

+1.250976562 

OCT 

5776377777777777 

DEC. VALUE  = 

-1.000000000 

OCT 

2004460200000000 

DEC. VALUE  = 

+9.507812500 

OCT 

1776462000000000 

DEC. VALUE  = 

+0.2988281250 

OCT 

2001500100000000 

DEC. VALUE  = 

+1.250488281 

OCT 

5776377707007777 

DEC. VALUE  = 

-1.000434756 

OCT 

2004460300000000 

DEC. VALUE  = 

+9.51 1718750 

OCT 

1776464300000000 

DEC. VALUE  = 

+0.3011474609 

OCT 

2001500000000000 

DEC. VALUE  = 

+1.250000000 

OCT 

5776377777077077 

DEC. VALUE  = 

-1.000006689 

OCT 

2004460400000000 

DEC. VALUE  = 

+9.515625000 

OCT 

1776465000000000 

DEC. VALUE  = 

+0.3017578125 

BSS 

1 

BSS 

1 

BSS 

1 

BSS 

1 

BSS 

2 

BSS 

1 

BSS 

2 

BSS 

1 

BSS 

2 

BSS 

1 

BSS 

3 

ENI 

3 

23 

RESET 

LMNT  INIT.  (SEE  K120) 

ENI 

0 

0 

LDA 

3 

ELI 

STA 

3 

ELI  INIT 

UP 

3 

/-I 

LDA 

0 

NT 

AJP 

1 

K20+5 

SLS 

0 

BEGIN 

HALT. 

START  RENDEVOUS  PROG. 

BSS 

3 

BSS 

1  17 

Bl=3 

66=65447  P=460 

BSS 

65 

B  1=6 

66=65566  P=460 

BSS 

1  16 

BUT 

66=65653  P=460 

BSS 

1 

STORE 

FOR  A/(1-WT) 

OCT 

2001552023631500 

OCT 

2002400000000000 

OCT 

2002600000000000 

OCT 

7777777777777777 

DEC 

lOOCO. 

OCT 

1754400000000000 

REM 

12 

RENDEVOUS  PROG.  STARTS  HERE 

ENI 

6 

SET  NUMBER  OF  ITERATIONS 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

LDA 

0 

V 

FDV 

0 

R 

V/R  = 

W 

STA 

0 

W 

STORE 

TARGET  ANGULAR  VEL- 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENA 

0 

35 

• 

STA 

ENI 

ENI 

0 

0 

0 

RUNKT 

0 

0 

SET  GILL  S.R.  FOR  29  DIFF.EQUS. 

ENI 

1 

0 

ENA 

0 

0 

CLEAR 

RUNKT4-17  AREA 

STA 

1 

RUNKT+17 

-  19  - 


L20 


L30 

RECUR 


LAMMU 


L50 


STA 

1 

RUNKT+20 

INI 

1 

2 

ISK 

1 

1  12 

SLJ 

0 

/-2 

LDA 

0 

UNITY 

LAMl (0)=! 

STA 

0 

RUNKT+20 

SET 

STA 

0 

RUNKT+37 

SET 

MU2 (0)=1 

STA 

0 

RUNKT+56 

SET 

PI3(0)=1 

STA 

0 

RUNKT+75 

SET 

RH04(0)=1 

STA 

0 

RUNKT+25 

SET 

PI]00T(0)=1 

LAC 

0 

UNITY 

STA 

0 

RUNKT+47 

SET 

LAM3D0T(0)= 

STA 

0 

RUNKT+66 

SET 

MU4D0T(0)=- 

LAC 

0 

TWO 

STA 

0 

RUNKT+44 

SET 

RH0200T(0)= 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENA 

0 

0 

STA 

0 

RUNKT+2 

SET 

TAU  =  0 

STA 

0 

RUNKT+4 

SET 

X (0)=0 

LDA 

0 

VI 

STA 

0 

RUNKT+3 

SET 

XD0T(0)=V1 

STA 

0 

RUNKT+12 

SET 

U( 0)=V1 

LOA 

0 

V2 

STA 

0 

RUNKT+6 

SET 

Y00T(0)=V2 

STA 

0 

RUNKT+ 15 

SET 

V( 0)=V2 

LDA 

0 

UNITY 

STA 

0 

RUNKT+7 

SET 

•< 

O 

11 

ENI 

0 

0 

SLJ 

0 

/+5 

BYPASS 

SIU 

6 

/+! 

SLS 

1 

/+! 

STOP  SWITCH 

ENI 

6 

0 

ENI 

1 

147 

LDA 

1 

SAVE 

STA 

1 

RUNKT  , 

UP 

1 

/-I  4 

ENI 

0 

0 

SLJ 

4 

EQUAL 

ZRO 

0 

0 

STA 

0 

RUNKT+1 

STORE  EQUALIZED 

SLJ 

4 

LAMMU 

RET. 

JUMP  TO  OBT 

SLJ 

0 

L70 

BYPASS  LAMMU 

ZRO 

0 

0 

SLJ 

0 

0 

ENTRANCE  TO  S.R 

SIL 

1 

L60+3 

SAVE 

:  Bi 

SIU 

3 

L60+4 

SAVE  83 

LDA 

0 

RUNKT+20 

LOAD  LAMl 

STA 

0 

ST0R3 

LAMl 

TO  ST0R3 

ENI 

0 

0 

ENI 

1 

0 

SET 

B  1=0 

ENI 

3 

0 

SET 

B3=0 

LDA 

1 

EL 

FMU 

3 

RUNKT+34 

EL(LAM2) 

FAD 

0 

ST0R3 

STA 

0 

ST0R3 

INI 

3 

14 

ENI 

0 

0 

• 

ISK 

1 

2 

SLJ 

0 

/-3 

STA 

0 

LAM 

STORE  LAMBDA 

FMU 

0 

LAM 

LAM 

SQUARED 

STA 

0 

STORl 

LAM 

SQUARE  TO  S 

LDA 

0 

RUNKT+23 

LOAD  MUl 

STA 

0 

ST0R3 

ENI 

3 

0 

LDA 

1 

EL 

FMU 

3 

RUNKT+37 

EL(MU2) 

-  20  - 


L60 


L70 


OERIV 


LUO 


L120 


FAD 

0 

ST0R3 

STA 

0 

STOR3 

INI 

3 

14 

6NI 

0 

0 

ISK 

1 

2 

SLJ 

0 

/-3 

STA 

0 

MU 

STORE  MU 

FMU 

0 

MU 

MU  SQUARED 

FAD 

0 

STORl 

LAM  SQ  +  MU  SQ 

STA 

0 

HYP2 

HYP2  =  LAM  SQ  +  MU  SQ 

SLJ 

4 

SQROOT 

ZRO 

0 

0 

STA 

0 

HYP 

HYP  =  SQ  ROOT  OF  HYP2 

FMU 

0 

HYP2 

STA 

0 

HYP3 

HYP3  =  HYP2  TIMES  HYP 

ENI 

1 

0 

RESTORE  B1 

ENI 

3 

0 

RESTORE  B3 

SLJ 

BSS 

0 

LAMMU 

3 

TO  EXIT  OF  S.R. 

LDA 

0 

A 

ENTER  FROM  L30+7 

FMU 

0 

LAM 

FDV 

0 

HYP 

STA 

0 

RUNKT+n 

SET  UOOTIO)  =  A(LAM)/HYP 

LDA 

0 

A 

FMU 

0 

MU 

FDV 

0 

HYP 

FSB 

0 

UNITY 

STA 

0 

RUNKT+14 

SET  VDOT(O) 

ENA 

0 

/+2 

CURRENT  ADDRESS+2  TO  A 

SAL 

0 

L220-4 

SET  EXIT  FROM  OERIV.  PROG. 

SLJ 

0 

L  150-3 

JUMP  TO  SET  AIJOOT(O) 

ENA 

0 

ALPHA+2 

SAL 

0 

L22C-4 

REPAIR  EXIT  FROM  OERIV. PROG. 

ENI 

0 

0 

SLJ 

0 

L220-3 

BYPASS  DERIV.  PROG. 

LAC 

0 

OMEGA 

FMU 

0 

RUNKT+2 

(-WT) 

FAD 

0 

UNITY 

1-WT 

STA 

0 

UNMINWT 

STORE  I 1-WT) 

LDA 

0 

A 

FDV 

0 

UNMINWT 

STA 

0 

ADOMINWT 

STORE  A/(1-WT) 

ENI 

0 

0 

LDA 

0 

RUNKT+12 

STA 

0 

RUNKT+3 

SET  XDOT  =  U 

LDA 

0 

RUNKT+15 

STA 

0 

RUNKT+6 

SET  YDOT  =  V 

LDA 

0 

RUNKT+4 

LOAD  X 

FMU 

0 

RUNKT+4 

X  SQUARED 

STA 

0 

ST0R2 

X  SO  TO  ST0R2 

LDA 

0 

RUNKT+7 

LOAD  Y 

FMU 

0 

RUNKT+7 

Y  SQUARED 

FAD 

0 

ST0R2 

X  SO  +  Y  SQ 

STA 

0 

STORl 

(X  SQ  +  Y  SQ)  TO  STORl 

SLJ 

4 

SOROOT 

FMU 

0 

STORl 

STA 

0 

R32 

STORE  RADIUS  CUBED 

FMU 

0 

STORl 

STA 

0 

R52 

STORE  RADIUS  FIFTH 

ENI 

1 

0 

SLJ 

4 

LAMMU 

NOTE  B1  SAVED  BY  LAMMU 

LAC 

\ 

RUNKT+4 

(-X)  TO  A 

FDV 

0 

R32 

(-X/R32)  TO  STORl 

STA 

0 

STORl 

LDA 

0 

ADOMINWT 

FMU 

1 

LAM 

FDV 

0 

HYP 

FAD 

0 

STORl 

STA 

1 

RUNKT+ 1 1 

SET  UOOT(T) 

INI 

1 

2 

21 


L130 


L140 


L150 


LJ60 


£NI 

0 

0 

ISK 

1 

5 

SLJ 

0 

L  12C-3 

£NI 

0 

0 

ENI 

3 

0 

LAC 

3 

RUNKT+26 

STA 

3 

RUNKT+  17 

LAC 

3 

RUNKT+3 1 

STA 

3 

RUNKT+22 

LOA 

0 

RUNKT+7 

FMU 

0 

RUNKT+7 

FSB 

0 

STOR2 

FSB 

0 

ST0R2 

FMU 

3 

RUNKT+20 

STA 

0 

STORl 

LAC 

0 

THREE 

FMU 

0 

RUNKT+4 

FMU 

0 

RUNKT+7 

FMU 

3 

RUNKT+23 

FAD 

0 

STORl 

FDV 

0 

R52 

STA 

3 

RUNKT+25 

LAC 

0 

RUNKT+7 

FMU 

0 

RUNKT+7 

FMU 

0 

TWO 

FAD 

0 

ST0R2 

FMU 

3 

RUNKT+23 

STA 

0 

STORl 

LAC 

0 

THREE 

FMU 

0 

RUNKT+4 

FMU 

0 

RUNKT+7 

FMU 

3 

RUNKT+20 

FAD 

0 

STORl 

FDV 

0 

R52 

STA 

3 

RUNKT+30 

INI 

3 

13 

ENI 

0 

0 

ISK 

3 

57 

SLJ 

0 

L  13C-4 

LAC 

0 

OMEGA 

FMU 

0 

RUNKT+2 

FAD 

0 

UNITY 

STA 

0 

UNMINWT 

LDA 

0 

A 

FDV 

0 

UNMINWT 

STA 

0 

AOOMINWT 

ENI 

0 

0 

LAC 

0 

MU 

FMU 

0 

RUNKT+20 

STA 

0 

STORl 

LDA 

0 

LAM 

FMU 

0 

RUNKT+23 

FAD 

0 

STORl 

STA 

0 

A1 

LAC 

0 

MU 

FMU 

0 

RUNKT+34 

STA 

0 

STORl 

LOA 

0 

LAM 

FMU 

0 

RUNKT+37 

FAD 

0 

STORl 

STA 

0 

A2 

FMU 

0 

A1 

FMU 

0 

ADOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+77 

LAC 

0 

MU 

FMU 

0 

RUNKT+50 

STA 

0 

STORl 

LDA 

0 

LAM 

FMU 

0 

RUNKT+53 

{-Pin 

SET  LAMIOOT  =  -PI  I 
(-RHO 1 ) 

SET  MUIOOT  =  -RHOl 

Y 

Y  SQUARED 

Y  SQ  -  X  SO 

Y  SQ  -  2{X  SQ) 

(Y  SQ  -  2  X  SQ)LAMl 

MINUS  THREE 
(-3X) 

(-3XY ) 

(-3XY  MUD 


SET  PI  IDOT 
(-Y) 

(-  Y  SQUARED) 

(-2  Y  SQUARED) 

X  SO  -  2  Y  SO 
(X  SQ  -  2  Y  SQ)MU1 

MINUS  THREE 
(-3X) 

(-3XY  ) 

I-3XY  LAMD 


SET  RH01D0T 


(-W).  ENTER  FROM  L70+5 
(-WT) 

1-WT 

LOAD  INITIAL  ACCEL. 

A/(  1-WT) 


{-MU  LAMl ) 


LAM  MUl 

A1  =  LAM  MUl  -  MU  LAMl 
{-MU  LAM2) 


LAM  MU2 

A2  =  LAM  MU2  -  MU  LAM2 
A1  TIMES  A2 
{A/  1-WT)A1A2 

SET  A12D0T 

{-MU  LAM3) 


LAM  HU3 


i 


22 


FAD 

0 

STORl 

STA 

0 

A3 

FMU 

0 

A1 

FMU 

0 

AOOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+102 

LAC 

0 

MU 

LI70 

FMU 

0 

RUNKT+64 

STA 

0 

STORl 

LDA 

0 

LAM 

FMU 

0 

RUNKT+67 

FAD 

0 

STORl 

STA 

0 

A4 

FMU 

0 

A1 

FMU 

0 

AOOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+ 105 

LDA 

0 

A2 

FMU 

0 

A2 

FMU 

0 

AOOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+1 10 

LDA 

0 

A2 

L200 

FMU 

0 

A3 

FMU 

0 

AOOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+ 1 13 

LDA 

0 

A2 

FMU 

0 

A4 

FMU 

0 

AOOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+ 1 16 

LDA 

0 

A3 

FMU 

0 

A3 

FMU 

0 

AOOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+121 

LDA 

0 

A3 

FMU 

0 

A4 

L210 

FMU 

0 

AOOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+1 24 

LDA 

0 

A4 

FMU 

0 

A4 

FMU 

0 

AOOMINWT 

FDV 

0 

HYP3 

STA 

0 

RUNKT+127 

ENI 

0 

0 

SLJ 

0 

ALPHA+2 

ENI 

0 

0 

ENI 

1 

147 

LDA 

1 

RUNKT 

STA 

\ 

SAVE 

UP 

1 

/-I 

SLJ 

4 

ALPHA 

L220 

ZRO 

0 

RUNKT 

ZRO 

0 

DERIV 

UP 

2 

/+! 

SLJ 

0 

/+3 

ENI 

0 

0 

SLJ 

4 

ALPHA+ 1 

ENI 

0 

0 

SLJ 

0 

/-2 

LDA 

0 

CAPXFIN 

FSB 

0 

RUNKT+4 

STA 

0 

C3 

FMU 

0 

C3 

STA 

0 

ST0R2 

LDA 

0 

CAPYFIN 

FSB 

0 

RUNKT+7 

A3  =  LAM  MU3  -  MU  LAM3 
A1  TIMES  A3 

SET  A  13D0T 
(-MU  LAM4) 

LAM  MU4 

A4  =  LAM  MU4  -  MU  LAM4 
A1  TIMES  A4 

SET  A14D0T 
A2  SQUARED 

SET  A22D0T 
A2  TIMES  A3 

SET  A2300T 
A2  TIMES  A4 

SET  A24D0T 
A3  SQUARED 

SET  A33D0T 
A3  TIMES  A4 

SET  A34D0T 
A4  SQUARED 

SET  A44D0T 

END  OF  DERIV  PROG. 

ENTER  FROM  L70+7 

STORE  INIT.  CONDITIONS 
FOR  NEXT  ITERATE 

SET  UP  GILL  ROUTINE 
PARAMETER  WORD 

JUMP  IF  82  NOT  =  ZERO 
82=0.  JUMP  TO  RMS  COMP. 
82  NOT  =  ZERO 
INTEGRATE  AGAIN 

JUMP  TO  TEST  82 

BEGIN  R.M.S.  COMPUTATION 

CAPXFIN  -  XFIN 

CAPYFIN  -  YFIN 


23 


L230 


L240 


L250 


L260 


L270 


STA 

FMl) 

FAD 

STA 

LOA 

FSB 

STA 

FMU 

FAD 

STA 

LDA 

FSB 

STA 

FMU 

FAD 

SLJ 

ZRO 

FOV 

STA 

ENI 

ENI 

LOA 

STA 

INI 

UP 

SLJ 

ZRO 

SLS 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

LDA 

STA 

LDA 

STA 

LDA 

STA 

LDA 

STA 

LOA 

STA 

STA 

LOA 

STA 

STA 

LDA 

STA 

LDA 

STA 

STA 

LDA 

STA 

ENI 

ENI 

SLJ 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 

ZRO 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

1 

3 

3 

1 

3 
1 

4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

3 

4 
0 
0 
0 
0 
0 
0 
0 
0 
0 


C4 

C4 

ST0R2 
ST0R2 
CAPUFIN 
RUNKT+ 12 
Cl 
C  1 

ST0R2 

STOR2 

CAPVFIN 

RUNKT+15 

C2 

C2 

STOR2 

SQROOT 

0 

TWO 

RMS 

17 

55 

RUNKT+20 
BBl  1 
-3 
/-I 

MATRIX 

0 

L250-3 

0 

0 

5 

4 

BBl  1 
0 
0 
0 

BBl  1+120 
RUMKT+100 
BBl  1  +  100 
RUNKT+103 
BBl  1+101 
RUNKT+106 
BBl  1+102 
RUNKT+1 1 1 
BBl  1+103 
RUNKT+1 14 
BBl  1  +  104 
BBl 1+106 
RUNKT+  1  17 
BBl  1+105 
BBl  1  +  11  1 
RUNKT+ 122 
BBl  1+107 
RUNKT+125 
BBl  1+1  10 
BBl 1+112 
RUNKT+ 130 
BBl  1+113 
0 
2 

MATRIX 

L270-2 

0 

0 

4 

4 

BBl  1  +  120 

h 

BBl  1+100 
3 


CAPUF IN  -  UFIN 


CAPVFIN  -  VFIN 


STORE  ROOT  MEAN  SQUARE 

SET  4  BY  4  MATRIX 
LAMIMUIPI IRHOI  =  B 
(1=1,2,3,41  AT  TAUFIN 

INVERT  THE  MATRIX 

SINGULAR  MATRIX  ALARM  HALT 


CODE  TO  INVERT  =  5 
L  =  4 

BBll  =  ADDRESS  OF  B  MATRIX 
N=0  NOT  USED  IN  INVERT 


M=0 

BBl 

A12 

SET 

A13 

SET 

A14 

SET 

A22 

SET 

A23 

SET 

SET 

A24 

SET 

SET 

A33 

SET 

A34 

SET 

SET 

A44 

SET 


NOT  USED  IN  INVERT 
1+120  =  ADDRESS  OF  INVERSE 


TO  SET  4  BY 
Dll  =  A12 

D12  =  A13 

D13  =  A14 

D21  =  A22 

022  =  A23 
D31  =  A32 

023  =  A24 
041  =  A42 

032  =  A33 

033  =  A34 
D42  =  A43 

043  =  A44 


3  D  MATRIX 


FOR  USE  6TH  WORD  HENCE 
MULT.  B  INVERSE  BY  0 
THIS  WORD  NOT  USED 


4  =  CODE  TO  MATRIC  MULTIPLY 
L  =  4  ROWS  IN  B  INVERSE 
BBl  1+120  =  ADDRESS  OF  B  INVERSE 
N  =  4  COLUMNS  IN  B  INVERSE 
BBll+100  =  ADDRESS  OF  0  MATRIX 
M  =  3  COLS.  IN  THE  MATRIX  PROD. 


24 


1300 


L310 


L320 


L330 


ZRO 

REM 

REM 

0 

881  1+  1 

LDA 

3 

88  1  1+12 

STA 

3 

842 

UP 

3 

/-I 

ENI 

3 

2 

LDA 

3 

BBl 1+7 

STA 

3 

B32 

UP 

3 

/-I 

ENI 

3 

2 

LDA 

3 

881  1  +  4 

STA 

3 

822 

UP 

3 

/-I 

LDA 

0 

RUNKT+ 1 1 

FSB 

0 

CAPUDFN 

STA 

0 

88  1  1 

LDA 

0 

RUNKT+  14 

FSB 

0 

CAPVDFN 

STA 

0 

821 

LDA 

0 

RUMKT+1 2 

FSB 

0 

CAPUFIN 

STA 

0 

831 

LDA 

0 

RUNKT+15 

FSB 

0 

CAPVFIN 

STA 

0 

841 

SLJ 

0 

L310 

STA 

0 

BBl  1 

ZRO 

0 

C4 

ENI 

0 

0 

SLJ 

4 

MATRIX 

SLS 

0 

L 31 0+1 

ZRO 

0 

0 

ZRO 

0 

0 

ZRO 

0 

5 

ZRO 

0 

4 

ZRO 

0 

881  1 

ZRO 

0 

0 

ZRO 

0 

0 

ZRO 

0 

0 

ZRO 

0 

8811+120 

ENI 

0 

0 

SLJ 

4 

MATRIX 

ZRO 

0 

L320-1 

ZRO 

0 

0 

ZRO 

0 

0 

ZRO 

0 

4 

ZRO 

0 

4 

ZRO 

0 

BBl 1+120 

ZRO 

0 

4 

ZRO 

0 

Cl 

ZRO 

0 

1 

ZRO 

0 

Cl  +  4 

ENI 

0 

0 

SLJ 

0 

/+2 

STA 

0 

BBl  1 

ZRO 

0 

Cl  +  7 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

UP 

6 

/  +  ! 

SLJ 

0 

FINAL 

SIU 

6 

/+! 

SLS 

2 

/+! 

ENI 

6 

0 

8611+ 1  =  ADDRESS  OF  MATRIX  PROD 
L273-306  SETS  UP  MATRIX  OF  SYSTEM  TO  B 
SOLVED  FOR  OELTAUFIN,  VAREL i  VAREM  AND  VA 


UDOTFIN 

UDOTFIN  -  CAPUDFN 

(OR 

PASS 

IF 

RENO.  POINT 

fri  7 

FIXED 

VDOTF IN 

VDOTFIN  -  CAPVDFN 

(OR 

PASS 

IF 

REND*  POINT 

K'V 

FIXED 

UFIN 

UFIN 

UFIN  -  CAPUFIN  (OR  PASS  IF  REND.  POINT  ^IXED) 


VFIN 

VEIN  -  CAPVFIN  (OR  PASS  IF  REND.  POINT  FIXED) 

BYPASS  PRINT  OUT 
PARAMETER  WORD  TO  PRINT 


INVERT  SYSTEM  MATRIX 
SINGULAR  MATRIX  ALARM  HALT 


5  *  CODE  TO  INVERT 
L  =  4 

8811  =  ADDRESS  OF  MATRIX 
N=0  NOT  USED 

M=0  NOT  USED 

8811+120  =  ADDRESS  OF  INVERSE 

MULT.  INVERSE  BY  COL.  MATRIX< 
THIS  WORD  NOT  USED 


4  =  CODE  TO  MATRIC  MULTIPLY 
L  =  4  ROWS  IN  INVERSE 
8811+120  =  ADDRESS  OF  INVERSE 
N  =  4  COLS.  IN  INVERSE 
Cl  =  ADDRESS  OF  COLUMN  MATRIX 
M=1  COL.  IN  THE  COL.  MATRIX 
Cl+4=  ADO.  OF  COL.  MATRIX  RESULT 

BYPASS  PRINT  OUT 
PARAMETER  WORD  TO  PRINT 


JUMP  IF  86  NOT  »  ZERO 

86  =  0.  JUMP  TO  FINAL  PROG. 

STOP  SWITCH 
RESET  86 


-  25  - 


ENI 

1 

3 

LOA 

1 

C  1  +  4 

STA 

] 

VARTAU 

STORE  DELTAUFIN, VAREL, VAREM.VAREN 

UP 

1 

/-I 

FAD 

0 

TAUFIN 

NEW  TAUFIN  FOR  NEXT  ITERATE 

STA 

0 

SAVE+143 

STORE  NEW  TAUFIN  (NOTE  L30+3 ) 

ENI 

0 

0 

OR  SLJ  0  L350+5  IF  REND.  POINT  FI 

L340 

FMU 

0 

W 

FAD 

0 

S 

STA 

0 

SPWT 

SLJ 

4 

TRIG+70 

STA 

0 

COSSPWT 

FMU 

0 

R 

STA 

0 

SAVE+136 

NEW  CAPYFIN 

LDA 

0 

SPWT 

SLJ 

4 

TRIG 

ZRO 

0 

0 

STA 

0 

SINSPWT 

FMU 

0 

R 

STA 

0 

SAVE+ 135 

NEW  CAPXFIN 

LDA 

0 

COSSPWT 

FMU 

0 

V 

STA 

0 

SAVE+ 137 

NEW  CAPUFIN 

L350 

LAC 

0 

SINSPWT 

FMU 

0 

V 

STA 

0 

SAVE+140 

NEW  CAPVFIN 

LAC 

0 

COSSPWT 

FMU 

0 

VSQDR 

STA 

0 

SAVE+142 

NEW  CAPVOFN 

LAC 

0 

SINSPWT 

FMU 

0 

VSQDR 

STA 

0 

SAVE+141 

NEW  CAPUDFN 

ENI 

0 

0 

ENI 

0 

0 

ENI 

] 

2 

LOA 

1 

Cl  +  5 

FAD 

1 

EL 

STA 

1 

EL 

STORE  NEW  EL,  EM  AND  EN 

UP 

1 

/-I 

L360 

SLJ 

4 

700C7 

PRINT 

ZRO 

0 

0 

SAU 

0 

MOOT 

DECIMAL  DUMP  WITH  PANEL 

ZRO 

0 

S 

i 

EXF 

0 

70 

CLEAR  ARITH.  ERRORS 

SLJ 

4 

70007 

PRINT 

STA 

0 

CAPXFIN 

DEC.  DUMP 

ZRO 

0 

LAM 

SLJ 

0 

RECUR 

TO  NEXT  ITERATION 

ZRO 

0 

0 

BSS 

3 

POLAR 

BSS 

200 

Bl=33  B6=66370  P=460 

ELI 

BSS 

1 

EMI 

ess 

1 

ENl 

BSS 

1 

T1 

BSS 

1 

EL2 

BSS 

1 

EM2 

BSS 

1 

EN2 

BSS 

1 

T2 

BSS 

1 

EL3 

BSS 

1 

EM3 

BSS 

1 

EN3 

BSS 

1 

T3 

BSS 

1 

EL4 

BSS 

1 

• 

EM4 

BSS 

1 

EN4 

BSS 

1 

T4 

BSS 

1 

ELS 

BSS 

1 

EM5 

BSS 

1 

EN5 

BSS 

1 

26 


T5 

BSS 

1 

BSS 

26 

STORE  FOR  MATRIX  PRODUCT 

FINAL 

ENI 

1 

147 

COMPUTES  POINTS  ON  FINAL  TRAJ. 

LIL 

2 

FIXIPLON 

LOA 

1 

SAVE 

STA 

1 

RUNKT 

RESET  INITIAL  CONDITIONS 

UP 

1 

/“I 

SLJ 

4 

ALPHA 

SET  UP  GILL  ROUTINE 

ZRO 

0 

RUNKT 

PARAMETER  WORD 

ZRO 

D 

DERIV 

ENI 

0 

0 

SLJ 

4 

/+4 

RET. JUMP  TO  DUMP  SUBROUTINE 

UP 

2 

/+1 

JUMP  IF  B2  NOT  =  ZERO 

SLJ 

0 

L70C+1 

B2=0.  JUMP  TO  REMAINDER  OF  FINAL  PROG 

L650 

ENI 

p 

0 

SLJ 

4 

ALPHA+1 

INTEGRATE  AGAIN 

ENI 

Q 

0 

4 

SLJ 

0 

/-3 

SLJ 

0 

0 

ENTRY  TO  DUMP  SUBROUTINE  ' 

LDA 

0 

RUNKT+2 

TAU 

ENI 

0 

0 

ENI 

0 

0 

STA 

0 

DUMP 

(DUMP)  =  TAU 

LDA 

0 

RUNKT+4 

X 

ENI 

0 

0 

ENI 

0 

0 

STA 

0 

DUMP+4 

(OUMP+4)  =  X 

LOA 

0 

RUNKT+7 

Y 

ENI 

0 

0 

ENI 

0 

0 

L660 

STA 

0 

OUMP+5 

(DUMP+5)  =  Y 

LDA 

0 

RUNKT+12 

U 

ENI 

0 

0 

ENI 

0 

0 

STA 

0 

DUMP+6 

(DUMP+6)  =  U 

LOA 

0 

RUNKT+15 

V 

ENI 

0 

0 

ENI 

0 

0 

STA 

0 

OUMP+7 

(OUMP+7)  =  V 

LDA 

0 

RUNKT+1 1 

UOOT 

ENI 

0 

0 

ENI 

0 

0 

STA 

0 

DUMP+10 

(DUMP+10)  =  UDOT 

LOA 

0 

RUNKT+ 14 

VDOT 

ENI 

0 

0 

ENI 

0 

0 

L670 

STA 

0 

ouMp+n 

(DUMP+ni  =  VDOT 

LDA 

0 

RUNKT+12 

U 

LDQ 

0 

RUNKT+15 

V 

SLJ 

4 

POLAR+130 

ENI 

0 

0 

ENI 

0 

0 

STQ 

0 

OUMP+12 

(DUMP+12)  =  ARCTAN(V/U) 

LDA 

0 

LAM 

LDQ 

0 

MU 

SLJ 

4 

POLAR+130 

ENI 

0 

0 

ENI 

0 

0 

STQ 

0 

DUMP+13 

(DUMP+13)  =  ARCTAN(MU/LAM) 

SLJ 

4 

70007 

PRINT 

STA 

0 

DUMP 

PARAMETER  WORD 

ZRO 

0 

DUMP+13 

L700 

SLJ 

0 

L650+2 

TO  EXIT  OF  DUMP  S.R. 

ZRO 

0 

0 

LOA 

0 

C3 

C3  =  CAPXFIN  -  XFIN 

ENI 

0 

0 

ENI 

0 

0 

STA 

0 

OUMP+14 

(OUMP+14)  =  CAPXFIN  -  XFIN 

LOA 

0 

C4 

C4  =  CAPYFIN  -  YFIN 

ENI 

0 

0 

27 


ENI 

0 

0 

STA 

0 

DUMP+15 

LDA 

0 

C  1 

ENi 

0 

0 

ENI 

0 

0 

STA 

0 

DUMP+16 

LOA 

0 

C2 

ENI 

0 

0 

L710 

ENI 

0 

0 

STA 

0 

DUMP+17 

SLJ 

4 

70007 

ZRO 

0 

0 

STA 

0 

DUMP+14 

ZRO 

0 

DUMP+17 

SLJ 

0 

L3570+5 

ZRO 

0 

0 

STORl 

BSS 

1 

ST0R2 

BSS 

1 

STOR3 

BSS 

1 

ST0R4 

BSS 

1 

L720 

DEC 

1. 

DELXl 

BSS 

1 

DELY  1 

BSS 

1 

DELUl 

BSS 

1 

OELVl 

BSS 

1 

DEC 

1  . 

DELX2 

BSS 

1 

DELY2 

BSS 

1 

DELU2 

BSS 

1 

DELV2 

BSS 

1 

DEC 

1  . 

DELX3 

BSS 

1 

OELY3 

BSS 

1 

DELU3 

BSS 

1 

DELV3 

BSS 

1 

DEC 

1. 

DELX4 

BSS 

1 

DELY4 

BSS 

1 

DELU4 

BSS 

1 

DELV4 

BSS 

1 

DEC 

1. 

DELX5 

BSS 

1 

DELY5 

BSS 

1 

0ELU5 

BSS 

1 

SAVE 

BSS 

150 

RUNKT 

BSS 

135 

CAPXFIN 

BSS 

1 

CAPYFIN 

BSS 

1 

CAPUFIN 

BSS 

1 

CAPVFIN 

BSS 

1 

CAPUOFN 

BSS 

1 

CAPVDFN 

BSS 

1 

TAUFIN 

BSS 

1 

LAM 

BSS 

1 

HYP 

BSS 

1 

HYP3 

BSS 

1 

MU 

BSS 

1 

HYP2 

BSS 

1 

UNITY 

DEC 

1  • 

W 

BSS 

1 

DELTAU 

BSS 

1 

SPWT 

BSS 

1 

COSSPWT 

BSS 

1 

SINSPWT 

BSS 

1 

R32 

BSS 

1 

RS2 

BSS 

1 

IMTSIGN 

FRACSIGN 

BSS 

1 

BSS 

1 

FLOiPLON 

BSS 

1 

FIXiPLON 

BSS 

1 

(DUMP+15)  =  CAPYFIN  - 
Cl  =  CAPUFIN  -  UFIN 


(DUMP-H6)  =  CAPUFIN  - 
C2  =  CAPVFIN  -  VFIN 

(OUMP+17)  =  CAPVFIN  - 
PRINT 

PARAMETER  WORD 

JUMP  TO  PRINT  B  FINAL 


ALSO  OELV5  STORE 
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YFIN 

UFIN 
VFIN  ^ 


THREDR2 

OCT 

2002417 

TREHALF 

DEC 

1.5 

DELTFLSG 

8SS 

1 

NT 

8SS 

1 

THETA 

flSS 

1 

VSTART 

DEC 

.65 

DELTHETA 

RSS 

1 

DELXFIN 

BSS 

1 

DELYFIN 

BSS 

1 

OELUFIN 

BSS 

1 

DELVFIN 

BSS 

1 

MOOT 

DEC 

.0025 

A 

BSS 

1 

QQ 

BSS 

1 

C 

DEC 

10000. 

EL 

BSS 

1 

EM 

BSS 

1 

£N 

BSS 

1 

S 

BSS 

1 

A1 

BSS 

1 

A2 

BSS 

1 

A3 

BSS 

1 

A4 

BSS 

1 

VARTAU 

BSS 

1 

VAREL 

BSS 

1 

VAREM 

BSS 

1 

VAREN 

8SS 

1 

ALPHA 

BSS 

134 

BSS 

14 

DUMP 

BSS 

20 

881  1 

BSS 

1 

812 

BSS 

1 

813 

BSS 

1 

814 

BSS 

1 

821 

BSS 

1 

822 

BSS 

1 

823 

BSS 

1 

824 

BSS 

1 

831 

BSS 

1 

832 

BSS 

1 

833 

BSS 

1 

834 

BSS 

1 

841 

BSS 

1 

842 

BSS 

1 

84  3 

BSS 

1 

844 

BSS 

1 

BSS 

120 

Cl 

BSS 

1 

C2 

BSS 

1 

C3 

BSS 

1 

C4 

BSS 

1 

BSS 

4 

LOG 

BSS 

100 

REM 

REM 

L2000 

BSS 

1440 

REM 

L3440 

ENA 

0 

0 

ENI 

1 

17 

STA 

1 

BBl  1 

UP 

1 

/ 

LDA 

0 

TAUFIN 

SLJ 

4 

TRIG 

STA 

0 

BBl  1 

LDA 

0 

TAUFIN 

SLJ 

4 

TRIG+70 

ZRO 

0 

0 

STA 

0 

B21 

LDA 

0 

TAUFIN 

FMU 

0 

R00T2 

6663160  THREE/ROOT2 


NON-DIMEN.  INITIAL  SPEED 


FRACTIONAL  MASS  LOSS  PER  SEC. 


ROCKET  NOZZLE  VEL.  IN  FT. /SEC, 


INIT.  TARGET  SECTOR  ANGLE 


Bl*30  86=67340  P=460  ’ 


5  TIMES  L  SQUARED  CELLS 


STORE  FOR  MATRIX  PRODUCT 
81=12  86=67700  P=460 

LOG  S.R.  ALARM  EXIT  TO  KllO 
(67703)  =  750  65110  000  00000 
USNPGS  GEN.  DUMP  81  =  32  P=2 1 

AID  TO  LINEARIZED  TRAJECTORY  PROG< 


CLEAR  88 II  MATRIX  AREA 


(8811 )  =  SINT- 


(821)  =  COST 


29 


SINHTR2 


L3450 


L3460 


L3470 


L3500 


L3510 


SLJ 

4 

EXP+71 

STA 

0 

832 

(832)  =  SINHTR2 

LDA 

0 

TAUFIN 

FMU 

0 

R00T2 

SLJ 

4 

EXP+107 

STA 

0 

B42 

(B42)  =  C0SHTR2 

LDA 

0 

SIN 

SINE  IMPULSE  INTEGRAL 

FDV 

0 

C 

STA 

0 

813 

LDA 

0 

COS 

COSINE  IMPULSE  INTEGRAL 

FDV 

0 

C 

STA 

0 

823 

LDA 

0 

SINH 

SINH  IMPULSE  INTEGRAL 

FDV 

0 

C 

STA 

0 

833 

LDA 

0 

COSH 

COSH  IMPULSE  INTEGRAL 

FDV 

0 

C 

STA 

0 

843 

LAC 

0 

CAPYFIN 

STA 

0 

8  14 

LAC 

0 

CAPVFIN 

STA 

0 

824 

SLJ 

0 

L357C+2 

JLNMP  TC  PATCH 

STA 

0 

834 

LDA 

0 

CAPUFIN 

STA 

0 

844 

LAC 

0 

DELXFIN 

STA 

0 

Cl 

Cl  =  ADD.  OF  MATRIX  COL. 

LAC 

0 

DELUFIN 

STA 

0 

C2 

LAC 

0 

DELYFIN 

FMU 

0 

R00T2 

STA 

0 

C3 

LAC 

0 

DELVFIN 

STA 

0 

C4 

ENA 

0 

/+2 

SAL 

0 

L320+4 

CREATE  LINEAR  EQU.  ROUTINE  EXIT 

SLJ 

0 

L310 

JUMP  TO  LINEAR  EQU.  ROUTINE 

ZRO 

0 

0 

ENA 

0 

L320+6 

SAL 

0 

L320+4 

REBUILD  RENDEVOUS  PROGRAM 

LDA 

0 

C  1  +  4 

DEL  VI 

FAD 

0 

V  1 

STA 

0 

ST0R1 

(ST0R1)  =  NEW  VI 

FMU 

0 

ST0R1 

STA 

0 

ST0R2 

(ST0R2)  =  NEW  VI  SQUARED 

LDA 

0 

C2+4 

DEL  V2 

FAD 

0 

V2 

STA 

0 

ST0R3 

(ST0R3)  =  NEW  V2 

FMU 

0 

ST0R3 

FAD 

0 

ST0R2 

NEW  VI  SQ.  +  NEW  V2  SQ. 

SLJ 

4 

SQROOT 

ZRO 

0 

0 

STA 

0 

VSTART 1 

NEW  VSTART 

LDA 

0 

ST0R1 

NEW  VI 

FDV 

0 

VSTART1 

COSINE  NEW  THETA 

SLJ 

4 

POLAR 

TO  INVERSE  COSINE  S.R. 

STA 

0 

THETA  1 

NEW  THETA 

LDA 

0 

C3+4 

DEL  C 

FAD 

0 

C 

NEW  C 

STA 

0 

EJECT 

(EJECT)  =  NEW  C 

LDA 

0 

C4+4 

DIFFERENTIAL  OF  B 

FAD 

0 

B 

NEW  B 

STA 

0 

SECTOR 

(SECTOR)  =  NEW  8 

SLJ 

4 

70007 

PRINT  RESULTS  DECIMAL 

STA 

0 

VSTART1 

PARAMETER  WORD 

ZRO 

0 

SECTOR 

SLJ 

4 

70007 

PRINT  RESULTS  OCTAL 

ZRO 

0 

0 

ZRO 

0 

VSTART  1 

PARAMETER  WORD 
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ZRO 

0 

SECTOR 

£NI 

r. 

0 

ENI 

0 

0 

LDA 

0 

VSTARTl 

RESUME  COMP.  OF  INIT.  CONDITIOI 

STA 

0 

VSTART 

SET  NEW  VSTART 

LDA 

0 

THETAl 

STA 

0 

THETA 

SET  NEW  THETA 

LDA 

0 

EJECT 

STA 

0 

C 

SET  NEW  C 

LDA 

0 

SECTOR 

STA 

0 

B 

SET  NEW  8 

SLJ 

4 

TRIG+70 

ZRO 

0 

0 

L3520 

STA 

0 

COSB 

FMU 

0 

R 

STA 

0 

CAPYFIN 

TARGET  COOR-  AT  RENDEVOUS 

LDA 

0 

B 

SLJ 

4 

TRIG 

ZRO 

0 

0 

STA 

0 

SINB 

FMU 

0 

R 

STA 

0 

CAPXFIN 

TARGET  COOR.  AT  RENDEVOUS 

LDA 

0 

COSB 

FMU 

0 

V 

STA 

0 

CAPUFIN 

TARGET  VEL.  COMP.  AT  RENDEVOUS 

LAC 

0 

SINB 

FMU 

0 

V 

STA 

0 

CAPVFIN 

TARGET  VEL.  COMP.  AT  RENDEVOUS 

LAC 

0 

COSB 

L3530 

FMU 

0 

VSODR 

STA 

0 

CAPVDFN 

TARGET  ACC.  COMP.  AT  RENDEVOUS 

LAC 

0 

SINB 

FMU 

0 

VSQDR 

STA 

0 

CAPUDFN 

TARGET  ACC.  COMP.  AT  RENDEVOUS 

LDA 

0 

THETA 

SLJ 

4 

TRIG+70 

COSINE  THETA 

ZRO 

0 

0 

FMU 

0 

VSTART 

STA 

0 

V  1 

INITIAL  HORIZONTAL  VEL. 

LDA 

0 

THETA 

SLJ 

4 

TRIG 

SINE  THETA 

FMU 

0 

VSTART 

STA 

0 

V2 

INITIAL  VERTICAL  VEL. 

LDA 

0 

C 

ROCKET  NOZZLE  VEL. 

FSB 

0 

TEN GRAND 

NOZZLE  VEL.  DIFFERENCE 

L3540 

AJP 

2 

/+! 

JUMP  IF  DIFF.  IS  POS. 

SCM 

0 

MASK 

ABSO.  VALUE  OF  DIFF. 

THS 

0 

UNITY 

EXIT  IF  ONE  GREATER  THAN  DIFF. 

SLJ 

0 

SECTOR+2 

NO.  JUMP  TO  ADJUST  MOOT  AND  C 

LDA 

0 

C 

YES.  LOAD  ROCKET  NOZZLE  VEL. 

FMU 

0 

MOOT 

FDV 

0 

GACCEL 

STA 

0 

A 

NON-DIMEN.  ACCELERATION 

LDO 

0 

ELS 

SLJ 

0 

KlOO+2 

RETURN  TO  LIN.  TRAJ.  PROG. 

VSTARTl 

BSS 

1 

TEMP.  STORE  FOR  NEW  VSTART 

THETAl 

BSS 

1 

TEMP.  STORE  FOR  NEW  THETA 

EJECT 

BSS 

1 

TEMP.  STORE  FOR  NEW  C 

SECTOR 

BSS 

1 

TEMP.  STORE  FOR  NEW  8 

SLS 

0 

SECTOR+1 

HALT  IF  1-WT  IS  NEGATIVE 

ZRO 

0 

0 

LAC 

0 

TAUFIN 

(-T) 

FMU 

0 

OMEGA 

(-WT)  . 

FAD 

0 

UNITY 

(  1-WT  ) 

AJP 

3 

SECTOR+1 

HALT  IF  1-WT  IS  NEGATIVE 

SLJ 

4 

LOG+66 

LOG(I-WT)  TO  BASE  E 

ZRO 

0 

0 

FMU 

0 

C 

FDV 

0 

TENGRAND 

SLJ 

4 

EXP 

31 


ZRO 

0 

0 

SCM 

0 

MASK 

FAD 

0 

UNITY 

L3560 

FOV 

0 

TAUFIN 

STA 

0 

OMEGA 

LOA 

0 

GACCEL 

FDV 

0 

REARTH 

SLJ 

4 

SOROOT 

ZRO 

0 

0 

FMU 

0 

OMEGA 

STA 

0 

MOOT 

LOA 

0 

TENGRANO 

STA 

0 

C 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

L3570 

ENI 

0 

0 

LOA 

0 

MOOT 

FMU 

0 

C 

FOV 

0 

GACCEL 

STA 

0 

A 

SLJ 

0 

K70 

LDA 

0 

CAPXF  IN 

FMU 

0 

RC0T2 

SLJ 

0 

L34d0+2 

ZRO 

0 

0 

SLJ 

4 

70007 

ZRO 

0 

0 

STA 

0 

UNITY 

ZRO 

0 

SPWT 

SLS 

0 

L3570+7 

ZRO 

0 

0 

END 

NEW  OMEGA 
STORE  NEW  OMEGA 


SHULER  ANGULAR  FREQ. 

NEW  MOOT 
STORE  NEW  MOOT 
LOAD  NEW  C 
STORE  NEW  C 


STORE  NEW  A 

JUMP  TO  LIN.  TRAJ.  PROG. 
PATCH  FCR  L3I46O+  I  L.  I  . 


ENTRY  FROM  L710+3 

PRINT  B  FINAL  DECIMAL 

PARAMETER  WORD 

FINAL  B  =  SPWT 

END  OF  NON-LINEAR  PROGRAM 
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genTA  7.U6  no.34 

Orbital  transfer  in  minimum  time. 
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